Nondestructive readout for a superconducting flux qubit 
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We present a new readout method for a superconducting flux qubit, based on the measurement 
of the Josephson inductance of a superconducting quantum interference device that is inductively 
coupled to the qubit. The intrinsic flux detection efficiency and back-action are suitable for a 
fast and nondestructive determination of the quantum state of the qubit, as needed for readout 
of multiple qubits in a quantum computer. We performed spectroscopy of a flux qubit and we 
measured relaxation times of the order of 80 ^ls. 

PACS numbers: 03.67.Lx, 85.25.Cp , 85.25.Dq 
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Suitably designed superconducting circuits, based on 
Josephson junctions, behave as quantum two level sys- 
tems. Due to scalability and flexibility in their design 
parameters, they are promising candidates for quantum 
bits (or qubits), which are the basic units in a quantum 
information processor 0]. In such circuitSjCoherent evo- 
lution for single qubits was observed 0,@,0,H,IE0' ^'^'i 
a conditional gate for two qubits was demonstrated j^. 

Flux qubits consist of a superconducting loop inter- 
rupted by one or more Josephson junctions. The basis 
states have opposite persistent current. The quantum 
state can be read out by measuring the generated mag- 
netic flux, using a DC superconducting quantum inter- 
ference device (DC-SQUID). The critical current of the 
SQUID depends on this flux and is usually measured by 
determining the maximum supercurrent, where the de- 
vice switches to a finite voltage. When the SQUID is 
in this voltage state it generates quasiparticles that later 
recombine with a burst of energy. It also radiates strong 
high frequency signals into the whole circuit with, in fu- 
ture, multiple qubits and readout devices. Inevitably, 
significant quantum information is destroyed apart from 
the consequence of reading out one qubit. In contrast, we 
now measured the SQUID critical current by determin- 
ing the Josephson inductance, without dissipation in the 
SQUID system. In this Letter, we discuss the intrinsic 
properties of this inductive readout and we present the 
first results of measurements on a flux qubit. 

In the experiments, we use a persistent current qubit 
(PCQ) 0. The PCQ is a flux qubit consisting of a 
small inductance superconducting loop interrupted by 
three Josephson junctions. Two of the three junctions 
are equal, characterized by the Josephson coupling en- 
ergy Ej and the charging energy Ec-, and the third one 
is smaller by a factor a. At low temperatures and appro- 
priate values of Ej > Ec and a, and with an external 
magnetic flux close to (2n H- l)<I'o/2 (n=integer) in 
the loop, the circuit behaves as a two level system. In 
the basis of two circulating current or flux states, the 
Hamiltonian reads 



iJ = i(e(T^ + Aoi) 



(1) 



in which ai{i = x,y,z) are the Pauli matrices and 
e = 2Ip{^q\j — {2n + l)<I>o/2). The minimum energy level 
splitting A and the maximum qubit persistent current 
Ip are determined by Ej, Ec and a. The two energy 
eigenstates (which are linear superpositions of the current 
states) each have an expectation value of the circulating 
current given by the derivative of the corresponding en- 
ergy eigenvalue with respect to (see Fig. QJi and b). 

A DC-SQUID magnetometer contains a loop inter- 
rupted by two Josephson junctions. The SQUID criti- 
cal current is Ic = 2Ico \ cos(7r$sQ/*I*o)|) where Ico is the 
critical current of a single junction and <&sq is the flux 
in the SQUID loop, which contains a contribution due 
to the coupled flux qubit. The voltage V and the super- 
current I of the SQUID are related to the difference in 
superconducting phase across the terminals of the DC- 
SQUID, 7, through the Josephson relations: V = 0o7 
and I ^ Ic sin(7), in which (po = <I>o/27r. For small varia- 
tions of the current around an average value Iq < Ic, the 
DC-SQUID can be described as an inductor, with the 
Josephson inductance Lj ~ cjjo/ \/lc ~ Iq- 

The Josephson inductance is measured by injecting an 
AC current Iac of frequency v, generating a voltage 
^ {Iac / Ic)^ot^- The injected current leads to a small 
magnetic flux in the qubit loop with signiflcant frequency 
components at 0, v and 2iy. To avoid qubit excitation v 
is taken well below the qubit level splitting. The qubit 
induced change in the inductance is small (~ 1% in our 
experiment). To enhance the measurement sensitivity, a 
resonant circuit is formed by adding a shunt capacitor to 
the SQUID. If the frequency v is close to the resonance 
frequency of this circuit, a small change in the Josephson 
inductance leads to a large change in the AC voltage, re- 
sulting in efficient state detection. Our method is similar 
to the RF - single electron transistor, used to measure 
charge qubits [3, but avoids the measurement of a dis- 
sipative circuit element. A readout scheme similar to 
ours was proposed for charge qubits 1 ^^'^ ^ detection 
scheme exploiting the non-linearity of ^(7) was shown 
in 11] . A PCQ was studied in by measuring the 
qubit loop susceptibility using a coupled linear resonant 
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circuit, with response time on microsecond time scales. 

Figure Ic shows the schematic diagram of the circuit 
comprising the qubit, DC-SQUID and microwave excita- 
tion Hne. The combined PCQ and SQUID system is fab- 
ricated on a single oxidized silicon substrate using elec- 
tron beam lithography and double-angle shadow evapo- 
ration of aluminum films. The design values for the qubit 
Josephson junctions are Ec/h = 4.6 GHz, Ej/Ec = 110 
and a = 0.75, which yields Ip — 660 nA. The areas of 
the PCQ and SQUID are, respectively, Aqh = 100 ^m^ 
and AsQ = 128 /im^. The chip containing the qubit and 
SQUID is attached to a printed circuit board on which 
the capacitor C and the resistors Ri, Rb and Rm are 
mounted. Capacitor C is positioned close to the SQUID 
to reduce the stray inductance Lg. Resistor Rb deter- 
mines the (AC and DC) driving current to the SQUID 
resonator. Applying a DC current and measuring the 
SQUID DC voltage via resistor Rm allows to determine 
the switching current. Applying an AC current to the 
SQUID resonator generates an AC voltage, which is cou- 
pled via the resistor Ri = 820 fl to the 50 fl input 
impedance of the low noise 0.5-1 GHz amplifier. With 
Rm,Rb » Ri » 50 n, the quality factor Q ~ 40 is 
mostly determined by the coupling resistor Ri. An ex- 
ternal superconducting coil generates the magnetic fiuxes 
^qb.ext, in the qubit loop and ^sQ,ext, in the SQUID loop. 
Qubit excitation is performed by a microwave magnetic 
fiux, using a line with mutual inductance Mmw=0.12 pH 
to the qubit. Appropriate noise filtering is employed in 
all signal lines. The printed circuit board is mounted 
inside a Cu box, anchored at the mixing chamber of a 
dilution refrigerator operating at T < 30 mK. 

We now discuss the requirements necessary for the 
readout of the flux qubit. The measurement process dis- 
turbs the qubit state, resulting in energy relaxation with 
a rate Tr and randomization of the phase of the state 
wavefunction or dephasing with a rate F^. For a mean- 
ingful measurement the relaxation time Tr — 1 /Tr needs 
to be much larger than the discrimination time Tdiscr, 
which is the time required to obtain enough information 
to infer the qubit state. It is defined as 

Td^scr = Sv/{Vi - Vo)^ (2) 

where Sv is the spectral density of the detector output 
voltage noise and Vq and Vi are the voltage values cor- 
responding to the qubit in the ground and in the excited 
state, respectively. We have calculated the decoherence 
rates for the case of moderate AC current driving. With 
the details due to be presented separately, the energy 
relaxation and the dephasing rate are given by 

F. = ^sin^e fc2(7o)[S'|(cooi + 27:iy) + S+{iJni - 27r,y)] 

(3) 

and 




FIG. 1: (a) Ground and excited state energy levels versus the 
magnetic flux in the PCQ loop, for Ej/Ec — 110, Ec/h = 4.6 
GHz and a = 0.75. (b) Expectation value of the qubit persis- 
tent current in the ground and excited state, (c) Diagram of 
the experimental setup. A waveguide is used to apply a small 
magnetic flux with amplitude Mmmlmw in the qubit loop. A 
coaxial line is used to apply the AC and DC bias current. The 
output AC voltage is amplifled at 4.2 K and room tempera- 
ture with gains Gi — 14 dB and G2 = 70 dB. The values of 
the components are Ls — 2.35 uH, C = 12 pF, Rb — 5.6 kfl, 
Rm = 11 kQ, and Ri = 820 Q.. 



Here tan 6* = A/e, woi = Ve^ -I- A'^/h and 5"+ is the 
Fourier transform of the symmetrized correlation func- 
tion of the SQUID phase operator 7. The coupling of 
the qubit to the SQUID fc(7o) is given by 



fc(7o) = MIpIctrcSQlO 



(5) 
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cos^^ fc2(7o)[5+(2^i.)]. 



(4) 



where M is the qubit-SQUID mutual inductance and 
Icirc,SQ is the SQUID circulating current, given by 
/co sin(7r$sQ/$o)- The amplitude of the SQUID phase 
oscillations 70 is proportional to Iac and thus the cou- 
pling between the detector and the qubit can be changed 
in-situ. The discrimination time, the relaxation time and 
the dephasing time = 1 /F^ are inversely proportional 
to 7q . This fact illustrates the tradeoff between obtaining 
information about the state of a quantum system and the 
disturbance introduced by the measurement process. In 
practice, the discrimination time is limited by the noise 
of the first amplification stage. An analysis of the circuit 
shown in Fig. ^ shows that efficient readout is possi- 
ble if an optimized cryogenic amplifier [l^ is used for 
readout. A second requirement for the detector is that 
there is no significant infiuence on the qubit during state 
preparation and manipulation, when no measurement is 
performed. This case has been discussed in ^3|> and for 
our circuit this influence is found to be negligible on the 
time scales necessary for qubit manipulation. 

An initial characterization of the SQUID and qubit 
was done by measuring the switching current of the DC- 
SQUID O. The dependence of the average SQUID 
switching current on the magnetic field, shown in Fig. 
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FIG. 2: (a) SQUID average switching current versus mag- 
netic field, (b) SQUID average switching current at step 2, 
with a linear background subtracted, (c) Frequency response 
of the readout circuit at a magnetic field corresponding to 
qubit step 2. (d) The response to magnetic flux changes for 
two operating frequencies, indicated in (c). 



had a period ABsq — 16.1 /iT. At a few equidistant 
positions separated by AB — 21 /iT the SQUID mod- 
ulation curve shows an increased slope. We attribute 
these features to the qubit generated flux, which changes 
sign when — {2n + l)$o/2. This is confirmed by 
the fact that AB/ABsq — Asg/Aqi,, indicating a pe- 
riod of 1 $0 in the PCQ loop, and the variation of the 
average switching current is in the direction of increasing 
flux (see Fig. ^). The average switching current of the 
SQUID around B = —63.7 /iT is plotted versus ^qb.ext in 
Fig.[2h- For the inductive measurements we concentrate 
on the qubit steps aX B = —0.6 //T and dX B ^ —63.7 
/iT, separated by 3<i>o (indicated as steps 1 and 2). 

Now, turning to the inductive method, we measured 
the amplitude of the SQUID AC voltage versus the fre- 
quency of the injected AC current. A resonance peak 
was found at a frequency I'o, dependent on ^SQ,ext- 
The width of this peak and the magnetic field depen- 
dence are consistent with the expected properties of the 
LC resonance for the circuit shown in Fig. From 
a fit of the dependence voi'^sQ^ext) using the relation 
M'^SQ,ext) = l/(27rv/(Ls + Lj{<^sQ,ext))G), we find 
Ls ~ 2.35 nH and the minimum value of the Joseph- 
son inductance Lj^nm — 0.31 nH. The obtained Lj^min 
is in agreement with the value of the maximum critical 
current of 1.0 /iA, as calculated using the value of the 
SQUID normal state resistance. 

Figure |2t shows the circuit resonance peak at a mag- 
netic field corresponding to the qubit step 2. Here, the 
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FIG. 3: (a) Spectroscopy measurements at — 17 GHz, 
for Tmw — 1 fis, Td = /is, and Tm — 2 fis, for Imw — (top 
curve) and for Imm ~ 7 /lA (bottom curve). The readout 
circuit has — 840 MHz and is operated at z/ = 855 MHz 
and 7o ~ 0.5. The data is plotted after averaging over flux 
intervals of 49 /i$o and 42 /i"I>o, respectively, and the two 
curves are oflset. (b) /miu vs flux in the PCQ loop (dots) and 
a linear fit through the origin (continuous line). 



circuit resonance frequency increases with increasing 
magnetic field. After fixing a SQUID driving frequency 
either smaller or larger than the resonance frequency, we 
vary the external magnetic field (see Fig. ^X)- Super- 
imposed on the smooth slope, the qubit step is clearly 
visible, similar to the data presented in Fig. ^p. 

As a next step we performed spectroscopy and mea- 
sured the qubit relaxation time. For this purpose we 
first apply a microwave burst of frequency fmw and du- 
ration Tmw to excite the qubit. After a delay time Td the 
SQUID detector is switched on by injecting the current 

1 AC a-ud the output voltage is measured for a time Tm- 
The time Tm is taken larger than the SQUID response 
time, Q/(2nh'o) ^ 8 ns. To improve statistics, we average 
over typically 1000 of these measurement sequences. 

Figure 13 shows the spectroscopy measurement results. 
In Fig. Ot, Vout is plotted against magnetic flux at step 
1, after substraction of the linear background. If mi- 
crowaves are applied, a peak and a dip are observed at 
positions symmetrically around the center of the step, 
with a height equal to almost half the step size. The 
peak and the dip correspond to the condition that the 
frequency of the applied MW burst matches the qubit 
energy level separation, which results in an incoherent 
mixture of the ground and excited states. The ampli- 
tude of the peaks and dips decreases with decreasing mi- 
crowave power (data not shown) and the position of the 
peaks depends on fmw, as shown in Fig. IJl^- We did not 
reliably identify peaks and dips for fmw < 8.5 GHz and 
therefore we could not determine A. A linear fit of the 
peak positions gives a value of the maximum persistent 
current Ip = 870 ± 30 nA; this agrees well with the value 
estimated from the amplitude of the qubit step, and is 
in reasonable agreement with the design value of 660 nA. 
We have performed the same measurements at qubit step 

2 and we obtained similar results. 

Both the spectroscopy peaks and dips display a weak 
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FIG. 4: (a) Spectroscopy peak at 17 GHz (plotted data is 
averaged over A"l> = 56.2 /i"l>o). (b) Spectroscopy peak area 
vs delay time and a fit with an exponential decay dependence. 



variation, periodic with the applied magnetic flux (see 
FigEK)- This period is the same in the peak and the dip 
and independent of microwave frequency and power. The 
position of the individual subpeaks changes linearly with 
/mtu) following the peaks and dips in which they are con- 
tained. In terms of qubit energies, the sub-peaks splitting 
A$5 corresponds to IpA^Jh = 880 ± 130 MHz,which 
matches the SQUID resonance frequency v within the ex- 
perimental uncertainty. This suggests that the observed 
substructure originates from transitions between the en- 
ergy levels of the coupled qubit-SQUID system. The co- 
herence time of the qubit, estimated from the width of the 
spectroscopy subpeaks, is T2 ~ 1 ns. This short decoher- 
ence time is not due to the measurement, as the SQUID 
AC driving is switched off during qubit excitation. The 
curves in Fig.|31i show a region, in the center of the step, 
characterized by a smaller slope. The width of this re- 
gion decreases with decreasing amplitude of the SQUID 
AC driving. We attribute it to the influence of the AC 
magnetic field generated by the measuring SQUID on the 
qubit. 

In order to determine the relaxation time, we have 
measured the spectroscopy peak using different time de- 
lays Td- Figure 0] shows the results obtained at qubit 
step 1 for a level separation of 17 GHz. The area of the 
Lorentzian curves fitting the peaks is plotted versus Td in 
Fig. ^p. A fit with an exponential decay gives the value 
of the relaxation time Ti — 77 ± 12 /is. The obtained 
value was the same for peak or dip and independent of 
the applied microwave power. 

From the measurements we estimate the parameters 
characterizing the readout efficiency, referred to the end 
of the amplification chain: ^/Sv — 4 fiV/\/lTz and 

- Vb ~ 800 nV, which imphes Tr/Tdiscr ^ "i-l- The 
measurement fidelity, defined as the probability to infer 
the qubit state correctly, based on the detector output 
voltage, is ~70%. Further improvement of the efficiency 
is possible by optimizing the SQUID circuitry. Using 
an on-chip circuit with small stray inductances and op- 
timized power transfer to the amplifier can reduce the 



measurement time by a factor of 100. Also, the present 
temperature T„ ~ 20 K of our amplifier, can be reduced 
to T„ ~ 3 K by using an optimized cryogenic amplifier or 
to r„ < 100 mK by using a SQUID amphfier with 
proportional decrease in the measurement time. 

The preliminary experiment of this paper, spec- 
troscopy measurements on a flux qubit, yielded the proof 
of principle for this method. The decoherence time was 
short, but feasible changes of the circuit can remedy this. 
The absence of the spurious back-action effects associ- 
ated with switching of Josephson junctions, the contin- 
uous nature of the detection and the in-situ tunability 
of the qubit-detector coupling allow for fundamental ex- 
perimental studies of the quantum measurement process 
and for correlation measurements on two qubits. 
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